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We study the groundstates of weakly interacting atomic Bose gases under conditions of rapid 
rotation. We present the results of large-scale exact diagonalisation studies on a periodic geometry 
(a torus) which allows studies of compressible states with broken translational symmetry. Focusing 
on filling factor = 2, we show a competition between the triangular vortex lattice, a quantum 
smectic state, and the incompressible k = A Read-Rezayi state. We discuss the corrections arising 
from finite size effects, and the likely behaviour for large system sizes. The Read-Rezayi state is 
stabilised by a moderate amount of additional dipolar interactions. 
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I. INTRODUCTION 

Under conditions of rapid rotation, ultra-cold gases of 
bosonic atoms are predicted to enter a very interesting 
regime of strongly-correlated many-particle physics.^ At 
sufficiently high vortex density, quantum fluctuations of 
the vortices can drive a quantum phase transition^i^i^ 
from the condensed vortex lattice phase into strongly- 
correlated liquid phases which are bosonic versions of 
conventional and unconventional fractional quantum Hall 
statesiii^i^'i The condition for quantum melting of the 
vortex lattice was first discussed in Refi, where the "fill- 
ing factor" , V, was shown to be the relevant control pa- 
rameter, and exact diagonalisation studies were used to 
show evidence of a triangular vortex lattice for > 6. 
(For a uniform system, v = N/N^ where TV is the num- 
ber of bosons and A^v the number of vortices.) 

While the limits of large^i^i^i^ and smaUp^i^i^ii filling 
factors are well understood, the nature of the ground- 
states in the regime of intermediate filling factors is 
less clear. This is a difficult regime to study theoret- 
ically, but is likely to show some of the most interest- 
ing new physics and will be the first to be accessed 
as experiments progress towards the strongly-correlated 
regime?'^"^*' . In Ref,2 evidence was presented for the ap- 
pearance of the sequence of Read-Rezayi (RR) statesii at 
u = k/2 with k integer, which are incompressible liquid 
states whose quasiparticle excitations obey non-abelian 
exchange statisticsJ^ Subsequent exact diagonalisation 
studies have shown that good convergence to the RR 
states at > 3/2 is not achieved with available system 
sizes. ^'li^ The role of the RR state ai v = 3/2 was signif- 
icantly clarified in Ref[l2 where it was shown that this 
state provides an excellent description of the groundstate 
when a small additional non-local interaction is intro- 
duced; similar effects have been reported for v — 2 and 
5/2 in the spherical geometryiii 

In this paper we study the groundstates of rotating 
bosons in the regime of intermediate filling factor using 
exact diagonalisation studies in the torus geometry up 
to system sizes much larger than those previously re- 



ported in this regime. We show that the groundstates 
involve a competition between compressible states with 
broken translational symmetry and the incompressible 
RR states. Our results identify a competing "smectic" 
phase, and point out the important role that finite size 
effects can have in stabilising this state both in numerical 
studies and in experimental systems. 



II. MODEL 

We consider a system of bosonic atoms of mass M con- 
fined to a harmonic trap with cylindrical symmetry, and 
with frequencies c^n and a;_L in the axial and transverse 
directions setting the lengthscales a||^j^ = ^J'h/ {Mlo\\ i). 
The atoms interact through contact interactions V{r) — 
g6^{r), with g — Anh^as/M chosen to reproduce the 
s-wave scattering length as. We study the regime of 
weak interactions, when the mean interaction energy is 
small compared to the trap energies, such that the par- 
ticles are restricted to the quasi-2D lowest Landau level 
(LLL) regime jiiSiii^ In this limit, the characteristic en- 
ergy scale is set by the s-wave Haldane pseudopoten- 
tial, Vq = 5/[(27r)'^/^a^a|[]. We have investigated the 
nature of the bulk groundstates {i.e. of systems con- 
taining large numbers of vorticesi^') using exact diago- 
nalisation studies on a periodic rectangular geometry (a 
torus) of sides a and 5, which contains iVv = ab/(TTa']_) 
vortices.— The torus lends itself naturally to the present 
study: it is commensurate with crystalline order so al- 
lows the study of states with broken translational invari- 
ance; the existence of a conserved momentum!^ reduces 
the sizes of the matrices in the diagonalisations and per- 
mits the identification of translational-symmetry-broken 
states through the appearance of quasi-degeneracies at 
wavevectors equal to the reciprocal lattice vectors (RLVs) 
of the crystalline order We express the momentum as 
a dimensionless wavevector K = (Kx,Ky) using units 
of 2TTh/a and 2'Kh/b for the x and y components, and 
report only positive Kx,Ky [states at {±Kx,iKy) are 
degenerate by symmetry]. 
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FIG. 1: Energy spectrum at z/ = 2 as a function of 
yJlCf^^^^K^, with torus aspect ratio chosen to be close to the 
value for which the groundstate energy is a local minimum, 
(a) iVv = 8, iV = 16, a/6 = 2 (b) A^v = 12, iV = 24, a/b = 4/3. 
The solid symbols mark the two lowest-energy states indi- 
cating the ordering wavevector of the smectic groundstate. 
The shaded symbols mark a band of excitations discussed in 



the text. [Note that the components 
units of 27r/a and 27r/6.] 



Ky are expressed in 



III. RESULTS 

We find that for filling factors v > 2 the spectrum is 
sensitive to the aspect ratio a/h of the torus.— The sensi- 
tivity to boundary conditions is a sign of a compressible 
groundstate, at least on the lengthscale defined by the 
finite system size in the calculation. Here, we shall focus 
on filling factor v = 2, which is convenient to study from 
a numerical point of view. (The qualitative features are 
reproduced at larger filling factors.) 



A. Smectic Groundstate 

In Fig. [1] we present the low-energy spectra &tv — 2 for 
systems containing A'v = 8 and 12 vortices, at aspect ra- 
tios (a/6 = 2 and 4/3) which are very close to the values 
for which the groundstate energy has a local minimum, 
In both cases, the groundstate [at K = (0, 0)] is quasi- 
degenerate with a single state at non-zero wavevector [at 



K = (4,0)]. This indicates a tendency to the formation 
of a smectic groundstate^i - i.e. a density wave state 
in which translational symmetry is broken in only one 
direction. In Figs[TJa) and (b) the smectic groundstates 
consist of four stripes parallel to the y-axis, with 4 and 
6 particles per stripe respectively. Although the smectic 
period is close to the lattice constant of a square lattice, 
a^'^ = y/TTa±, there is no indication of ordering in a vor- 
tex lattice (there are no quasi-degenerate states at the 
requisite RLVs). 

The low-energy excitation spectra above these quasi- 
degenerate states are consistent with the existence of a 
smectic groundstate. For A^v = 8, Fig. [TJa), the ground- 
state has four stripes of four particles each. There exist 
low-energy "particle-hole" excitations in which a particle 
is moved from one stripe to another (to give occupations 
of 4,4,3,5). These excitations account for the band of 
states shown as shaded symbols in Fig. [1] (a) . This band 
is narrow showing that exchange interactions between the 
stripes are small. For larger systems (longer stripes) one 
expects the energy of this band to fall and exchange inter- 
actions to increase, so these particle-hole excitations will 
contribute increasingly to the groundstate. It is possible 
that this will lead to a phase locking of the stripes, caus- 
ing the formation of a vortex lattice: a set of states at 
the RLVs of the lattice would then emerge from this band 
and become quasi-degenerate with the groundstate; the 
remaining states would contribute to the phonon mode of 
the lattice. The exact diagonalisation studies show that, 
at least up to = 12, a vortex lattice does not form. 
Rather, the system shows behaviour expected of a smec- 
tic groundstate: the analogous low-energy band [shaded 
symbols in Fig. Wlh)] broadens to give low-energy states 
at K = (0,3) and (0,6) corresponding to the competing 
re-oriented smectic state (with three stripes of eight par- 
ticles parallel to the x-axis); the remaining low-energy 
states in this band form the expected gapless shear mode 
of this smectic state. Thus, the exact diagonalisation re- 
sults up to iVv = 12 are consistent with a smectic ground- 
state at t/ = 2.— 



B. Triangular Vortex Lattice 

Evidence for a competing vortex lattice state aX v — 2 
can be obtained by choosing an aspect ratio that is com- 
mensurate with the triangular vortex lattice (away from 
the local minimum in groundstate energy). There then 
develop approximate quasi-degeneracies of the ground- 
state with states at the RLVs of the triangular lattice, 
which emerge from the band of low-energy states dis- 
cussed above. The spectra are shown in Figl2] for sys- 
tems with Av — 8, 12. In all systems we have studied the 
states at the RLVs of the triangular lattice do not become 
well separated from other states in the low-energy band. 
Moreover, since this aspect ratio does not correspond to 
a local minimum of the energy, there is poor evidence 
that the groundstate is a triangular vortex lattice. 
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FIG. 2; Energy spectrum at z/ — 2 as a function of 
Ki + Ky , with torus aspect ratio chosen to be consistent 
with a triangular vortex lattice, (a) A'v = 8, = 16, a/6 — 
(b) TVv = 12,iV = 24, a/6 = 2/V3. The solid symbols 
mark the wavevectors of a triangular vortex lattice. These 
states are not well separated from higher energy states, show- 
ing a weak tendency to vortex lattice ordering. [Note that the 
components Kx,Ky axe expressed in units of 27r/a and 27r/6.] 



C. Finite-size Effects 

As we now discuss, finite size corrections may be acting 
to favour the smectic state as compared to the triangular 
vortex lattice. We illustrate the relevance of finite size 
corrections by making simple mean-field ansatze for the 
triangular vortex lattice and smectic states. 

The triangular vortex lattice is the fully condensed 
state in the 2D LLL that minimises the mean interac- 
tion energy. The mean energy of this state on a periodic 
geometry is 



(1) 



where /3a = 1.1596^?^ and the correction at finite iVv 
arises from the lack of self-interaction of the particles [the 
mean interaction energy is proportional to N(N — 1)]. 

As a simple description of the smectic state, we con- 
sider the Fock state, defined on a long cylinder with cir- 
cumference Ly^ in which iVg bosons occupy each of a set 
of Landau-gauge orbitals ("stripes") that are extended 



around the circumference of the cylinder and that are 
equally spaced along the length of the cylinder by Ax. 
(This state is similar to the Tao-Thouless statei^) The 
smectic state describes a fragmented condensate, with Ns 
particles in each stripe and with no fixed phase relation- 
ship between the stripes. (In the LLL, the establishment 
of a fixed phase relationship between the stripes would 
establish fixed positions of the vortices, i.e. form a vor- 
tex lattice.) For a large circumference, Ly oo, the 
mean interaction energy is minimised, at fixed filling fac- 
tor I' — NsTTa\/{LyAx), for a period Ax — 1.73a^. This 
is close to the period of the smectic state appearing in 
numerics, of approximately a^'^ = 1.77a±. At this period 
the energy per particle of the smectic state is 



^smectic (Ax/a^'i) 

= lyfJsVo Vq 



N 



(2) 



where f3s = 1.17184, and A^v.s = Ng/i' is the number 
of vortices associated with each stripe. The finite-size 
correction arises from the lack of self-interaction between 
the Ns particles in each stripe. 

Since (3s > /9a, the smectic state has higher energy 
than the triangular lattice in the thermodynamic limit 
(iVv, iVv,s oo). This is a consequence of the general re- 
sult that repulsive interactions disfavour fragmentation of 
a condensate owing to the loss of exchange interactions.— 
It is remarkable, however, just how small the energy dif- 
ference between these states is [(/3s — Pa) /Pa !%]■ 
For a system on a long cylinder with finite circumference 
(iVv oo with iVv,s fixed) the smectic state has lower en- 
ergy than the triangular vortex lattice for iVv,s ^ 80/j^. 
The stabilisation of the smectic state can be viewed as 
a form of Mott transition: by fixing a definite particle 
number on each stripe, there is a reduction in Hartree 
energy at the expense of an increase in exchange energy. 
At J' = 2 these finite size effects cause the trial smectic 
state to have lower energy than the triangular vortex lat- 
tice up to systems of 40 vortices per stripe. This is far 
beyond the system sizes that can be studied in exact diag- 
onalisations. It is also far beyond the maximum number 
of vortices crossing the clouds in existing experiments^>i^ 
(on the order of ten), showing that finite size effects can 
play a significant role in determining the groundstates at 
intermediate filling factors in experimental geometries. 

The simple mean-field analysis presented here is not 
intended to provide quantitative corrections to the exact 
diagonalisation results. This analysis neglects quantum 
fluctuations which reduce the mean energies of the smec- 
tic and vortex lattice states even in the thermodynamic 
limit. Nevertheless, there are indications that finite size 
corrections do play an important role in the energetics of 
the smectic state in the numerics: the smectic state has 
a lower energy the shorter is the length of the stripes. 
For example, for = 12, a/6 = 4/3 [Fig. [TJb)] the two 
high symmetry configurations correspond to configura- 
tions with either four short stripes with iVv.s — 3 parallel 
to the X-axis, or three long stripes with 7Vv,s = 4 paral- 
lel to the y-axis; they have energies [obtained from the 
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respective ordering wavevectors, K = (4, 0) and (0, 3)] 
of E/N = 1.786Vb and E/N = 1.796^0- Similarly, 
from results on A'v — 12, a/6 — 3 one can obtain an 
estimate of the energy for short stripes with A^v,s = 2 
{E/N ~ 1.739Vb). There is a monotonic increase in the 
energy of the smectic state as iVv,s increases. Owing to 
the small number of system sizes available, it is difficult 
to extrapolate to the thermodynamic limit for either the 
smectic or vortex lattice state. However, the sizes of 
the energy changes of the smectic are comparable to the 
separation in energy to the triangular vortex lattice, for 
which estimates of the energy per particle can be as low 
as E/N = l.lSOVbi^ Thus, while the smectic state is the 
groundstate in systems of up to A^v = 12, it is possible 
that finite size effects are playing a role in stabilising the 
smectic state. 

In the studies of Ref 2 in which a transition to a vor- 
tex lattice was seen in systems of A^v = 6,8 vortices, 
the aspect ratio used was close to the minimum of the 
groundstate energy. In view of the above analysis, we 
can now say that in those studies^ the estimate of the 
critical filling factor ~ 6 arose from a competition be- 
tween the triangular vortex lattice and the smectic state. 
The finite-size corrections to the energy of the state com- 
peting with the triangular vortex lattice were not recog- 
nised in Ref,'2'. From the present analysis, we conclude 
that the analysis of Ref,2 underestimated the stability 
of the vortex lattice in the thermodynamic limit. The 
available numerical results are consistent with the vortex 
lattice being more stable than the smectic state down to 
as small a filling factor as — 2. 



D. Read-Rezayi State 

An interesting phase that competes with the smectic 
and the triangular vortex lattice is the k — 4 Read-Rezayi 
statejii As discussed above, for contact interactions the 
groundstate at = 2 appears to be compressible on the 
scale of the systems studied. However, as for ly — 3/2^^ 
by introducing interactions with non-zero range, we find 
that the spectrum becomes characteristic of the incom- 
pressible RR state (showing clear quasi-degeneracies at 
the expected wavevectors, and a gap to other excita- 
tions). The case of dipolar interactions^! is of partic- 
ular experimental relevance, as these can be important 
for atomic species with large magnetic moments»22i Fig. [3] 
shows the overlap of the exact groundstate with the RR 
state as a function of the additional dipolar interaction. 
In a narrow range around V2/V0 — 0.3 the groundstate 
is well described by the fc 4 RR state. For V2/V0 > 0.4 
the groundstate is well described by broken-symmetry 
states with stripe and "bubble crystal" ordering. These 
states have periods that are larger than a*"^, and are much 
more stable to quantum fluctuations^^ than are the com- 
pressible states discussed above. An enhancement of the 
overlap with the fc = 4 RR state in related calculations 
on the spherical geometry for a model short-range inter- 
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FIG. 3: (Color online) Variation of the overlap (squared) of 
the groundstate for contact and dipolar interactions with the 
fc = 4 RR state, for (a) iVv = 8, iV = 16, a/6 = 1 and (b) A^^ = 
10, A*' = 20, a/b = 1, for a±/a^^ <^ 1. We parameterise the 
relative size of dipolar and contact interactions by the ratio 
of the first two Haldane pseudopotentials V2/V0, as defined 
in Refs [l^i27l . The k = 4 RR state has five degenerate states 
on the torus located at the zone centre and cornersj^iii for 
K = (0, 0) the RR groundstate is a doublet and we plot the 
overlap of the lowest two states within this subspace, so the 
maximum overlap squared is 2. [In (b), the discontinuity of 
the overlap for K — (0,0) close to V2/V0 = 0.13 is due to a 
level crossing with a state of different point-group symmetry.] 
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IV. SUMMARY 

We have shown that exact diagonalisation studies of ro- 
tating Bose gases with contact interactions indicate that 
the groundstate at = 2 is a quantum smectic state. 
We argued that finite size corrections may lead to an en- 
hanced stability of the smectic state in the system sizes 
amenable to exact diagonalisations, and the triangular 
vortex lattice may have lower energy in the thermody- 
namic limit. In the regime of intermediate filling fac- 
tors, ly > 2, the correlation length of the groundstate 
for contact interactions remains larger than the maxi- 
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mum system sizes that are available in exact diagonal- 
isations. Comparisons with the RR states ei% v — k/2 
quickly become severely limited, as these states involve 
clustering of groups of k particles {d.t v = 2 exact diag- 
onalisations currently allow studies of at most six clus- 
ters). The torus is likely to favour broken translational- 
symmetry states of short period, so it is possible that 
the true groundstates at intermediate filling factors are 
quantum states of more complex symmetry - including 
quantum nematic phases^ or incompressible liquids with 



small gaps. Clearly this is a regime in which there are 
strong quantum fluctuations of vortices, and where novel 
strongly correlated groundstates may form. It will be 
very interesting to explore this regime experimentally. 
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(E.H.R.). 
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